On the frequency of gravitational waves 
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We show that there are physically relevant situations where gravitational waves do not inherit 
the frequency spectrum of their source but its wavenumber spectrum. 
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I. INTRODUCTION 



Let us consider a source of (weak) gravitational waves 
(GWs), an anisotropic stress IIy-(x,t). The induced 
GWs can be calculated via the linearized Einstein equa- 
tions. In the case of an isolated, non-relativistic, far 
away source, one can derive the quadrupole formula (see 
e-g- Q), r = |x|,c = 1 



2G 

hij(x,t) = —Qij(t - r) 



(1) 



where Qij denotes the quadrupole of the source, and we 
are considering the perturbed metric = r\^ v + h^ v . 
The wave has the same time dependence as the source. 
If we have an isolated source which has a harmonic time 
dependence with frequency u) s , IIjj(x, t) = Uij(x)e~ luJst , 
the wave zone approximation gives, far away from the 
source, 



hij (x, t) 



AGe - iu)s (t-r) 



d 3 x' IIy(x')e" 



(2) 



Again, hij has inherited the frequency of the source: it 
is a spherical wave whose amplitude in direction x is de- 
termined by Iljj(k) with k = cj s x. As we will show in 
this brief report, this simple and well known fact from 
linearized general relativity and, equivalently, electro- 
dynamics, has led to some errors when applied to GW 
sources of cosmological origin. 

As an example we consider a first order phase tran- 
sition in the early universe. This can lead to a pe- 
riod of turbulent motion in the broken phase fluid, giv- 
ing rise to a GW signal which is in principle observ- 
able by theplanned space interferometers LISA or BBO 
@, d, EI IE IE S Hi • One can describe this phase of tur- 
bulence in the fluid as a superposition of turbulent ed- 
dies: eddies of characteristic size I rotate with frequency 
lui ~ v 1 '£ ^ l/£, since v < 1. The wavenumber spectrum 
of the eddies peaks at fcj, ~ 1/L, where L is the stirring 
scale, and has the usual Kolmogorov shape for k > 1/L 
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(£ < L). The above considerations about GW genera- 
tion prompted several workers in the field 0, [1, 0, S @| 
to infer that the induced GWs will directly inherit the 
frequency of the eddies. Especially, even if the turbulent 
spectrum in A:-space peaks at wave number fc^ ~ 1/L, 
Refs. [2, S H B, H conclude that the GW energy den- 
sity spectrum peaks instead at the largest eddy frequency 
k = LdL 7^ kt,' The argument to reach this conclusion 
is simple and seems, at first sight, convincing: one may 
consider each eddy as a source in the wave zone which os- 
cillates at a given frequency. The GWs it produces have 
therefore the same frequency. Adding the signal from all 
the eddies incoherently, one obtains a GW energy density 
spectrum which peaks at the same frequency as the fre- 
quency spectrum of the eddies, namely at u = ujl- Since 
GWs obey the dispersion relation k = |cj|, the /c-space 
distribution of the GWs then cannot agree with the k- 
space distribution of the source. In this brief report we 
show why this simple argument is not correct. In the 
above case, the GW energy spectrum actually peaks at 
lu = ki, and the A:-space distribution (but not the fre- 
quency distribution) of the source and the GWs agree 
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Thus, the result @ has to be used with care in the 
case of a cosmological GW source which is active only 
for a finite amount of time and is spread over the en- 
tire universe. In Section IIIII we discuss several cases of 
cosmological sources to which this consideration applies. 
This is relevant especially in view of future searches for 
a stochastic GW background which might be able to de- 
tect some of the features predicted for GW spectra of 
primordial origin (in particular, in the case of turbulence 
after a phase transition, the peak frequency) Q . 

The rest of the paper is organized as follows. In Sec- 
tion II we analyze an academic example of a source which 
can be represented as a Gaussian wave packet both, in 
/c-space and in frequency space. We show when the wave 
will inherit the wave number of the source and when its 
frequency. We then discuss the case of a generic source. 
In Section III we consider a typical cosmological source: 
stochastic, statistically homogeneous and isotropic, and 
short lived. There we also explain the subtleties which 
are responsible for the above mentioned misconception 
in the case of turbulence, and we discuss other examples 
of GW sources for which our findings are relevant. The 
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final section contains our conclusions. 



II. THE GAUSSIAN PACKET 

Let us consider a scalar field /i(x, t) fulfilling the usual 
wave equation with source II (for the sake of simplicity, 
we omit the tensorial structure which is irrelevant for the 
problem we want to address here): 



□ft(x,t) = -167rGn(x,t) 



(3) 



Here □ = — d? + A. In terms of the Fourier transform of 



the source, 

fi(k,w)= Jdtj d 3 xe i(k x - wt) n(x,i), 
the retarded solution of Eq. <(3j) is (fc = |k|) 



(4) 



h{x,t) 



%G 

7T 2 



d A k- 



k 



n(k,-fc)e" 



ikt 



n(k,fc). 



i k / 



(5) 

To illustrate the possible behavior of this solution, we 
now consider a source whose Fourier transform has a 
Gaussian profile both in wavenumber and frequency, 1 



II(k, u) — Hp exp 



r (k~k s f- 




[ 2A 2 _ 


exp 



2A 2 



(6) 

where k s and u> s represent the characteristic frequency 
and wavenumber of this source (fc s is positive, while cu s 
can be either positive or negative). Inserting ^ in solu- 
tion one immediately obtains 



/i(x,t) 



-ikt 



exp 



AG n 



(k 



dk sin(fer) exp 



(k 



Y2l 



2A 2 



e ikt exp 



(fc - Usf 

2A 2 , 



where r = |x|. The resulting wave is a superposition of 
plane waves with frequency ±k. To estimate which is 
the dominant frequency of this superposition, i.e. the 
frequency at which the amplitude peaks, one has to look 
at the maximum of the double Gaussian in l[7]). In par- 
ticular for a source which has A w <C A&, we expect the 
spectrum to be peaked at k ~ ui s : we recover result (J2J, 
namely a spherical wave that oscillates with the char- 
acteristic frequency of the source. In the opposite case 
A w 3> Afe, the spectrum is peaked at k ~ k s , and the 
wave inherits the wavenumber of the source. 

This simple academic example provides us insight 
about the behavior of more realistic sources. Note that 



the spacetime distribution of the above wave-packet is 
again a Gaussian packet with spatial extension A x — 
1/Afc and temporal spread At = 1/A U . A typical as- 
trophysical source of GWs is long lived and has a very 
distinct frequency (supernovae, binary systems). On the 
other hand it is well localized in space, so we can expect 
the ordinary A w <C A& case. Note however, that the 
signal is suppressed if k s is very far from ±w s . For a 
typical binary we have w s ~ vk s , where v is the velocity 
of the system and it is well known that fast binaries emit 
a stronger GW signal than slow ones. 

On the other hand, consider a cosmological source 
which is spread over all of space in a statistically ho- 
mogeneous and isotropic way, to respect the symmetries 
of the Friedmann universe. This source is very extended 
but it may be active only over a short period of time. In 
a Friedmann universe, sources can be switched on and 
off everywhere at approximately the same time, e.g. at 
a given temperature. This kind of sources may then fall 
in the A w 3> A^ range. In the extremal case of a source 
which is a delta function in time, A w — > oo, the fre- 
quency spectrum of the generated radiation is entirely 
determined by the structure of the source in fc-space. 

The goal of the next section is to present an example 
of such a source. As a preparation, let us first consider 
a generic source II(x, t) which is active from some initial 
time ii n to a final time ifi n . Using the retarded Green's 
function in fc-space 



G R (k,t) = 47r9(i) 



sin(fci) 



(8) 



one finds the solution to Eq. ([3]) at time t > ts. n 



h(k, t) = A Sm{k( \ tfin)) + B cos(k(t - t fin )) , (9) 



with 



A(k) = 16ttG J ^ dt' cos(fc(t fin - t')) n (k, t 1 ) 

J tm k 



(10) 



If the source has some typical frequency uj s satisfying 
(tfin ~ ^in) w s 3* 1) the integrals A and B approximate 
delta functions of k — uj s and k + uj s . However, frequencies 
of the order of (tg n — tin) or smaller are not especially 
amplified and only a detailed calculation can tell where 
the resulting GW spectrum peaks. 



III. A STOCHASTIC, HOMOGENEOUS AND 
ISOTROPIC, SHORT LIVED SOURCE 



1 The source (0 is not real in physical space; to make it real 
one should symmetrize it in ui s . But we can as well consider a 
complex source whose real part represents the physical source. 



We now analyze the case of a typical GW source of 
primordial origin: a stochastic field, statistically homoge- 
neous and isotropic. Let us consider the two point func- 
tion of the spatial Fourier transform of the anisotropic 
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stress tensor 



(n iJ -(k,t)n^(k',t')} = 5(k-k')M iJ - H (fc,t,t'), (ii) 

and analogously for the induced GW 

(h %3 Ck,t)h* kl (k',t)) =5(k--k')H ijkl {k,t), (12) 

where we write (fc, t) for H^m (fc, t, t) . For simplicity 
we consider the wave equation ^ in flat space-time, and 
ignore the tensor structure of H and M . From the above 
definitions and from Eq. (|3|), one finds that in Fourier 
space the relation between the spectral power of the en- 
ergy momentum tensor and that of the GW is 



H{k,LU,u}') = (16irGY 



M(k,LO,Lo') 



(lo 2 - k 2 ){uj' 2 - k' z ) 



(13) 



Fourier transforming over lo and lo' and taking into ac- 
count the correct boundary and initial conditions H — > 
for t — oo, one obtains 



H(k,t) 



k' 2 



Re[M(k,k,k) - e- 2ikt M(k,k,~k)] , 

(14) 

where we have used that M(k,Lo, lo') = M*(k, —ui, —lo'). 
This result can also be obtained via convolution of 
M(k,t,t') with the Green's functions Gn(k,t)Gii(k,t') 
of Eq. ([5]) . To interpret Eq. (fT4"]) let us consider a source 
which is active from some initial time ti n to a final time 
ifi n and is oscillating with frequency lo s , 



M (M ,0 = {^ (t - n 



if 
if 



t e [tt 

t £ [tr 



(15) 



Inserting the Fourier transform of this somewhat simplis- 
tic source in the above Eq. (Tj"4"|) one finds 



H(k,t) = (WttG) 



2 F(k) sin 2 [(fc - Lo s )(t un 



iin)/2] 



k 2 



(k 



LO s 



(16) 

plus terms which average to zero over an oscillation pe- 
riod. The square root of the last factor in this expression 
becomes ir5(k — lo s ) when io s (t nn — t- m ) — > oo. Hence, if 
w s (ifin — tin) ^S> 1, the resulting GW spectrum H(k,t) is 
peaked at lo s , the typical frequency of the source. How- 
ever, if Lo s (t nn — ti n ) ~ C(l) ; the last term does not 
significantly influence the spectrum, which can then be 
peaked there where F(k) is. The typical frequency of 
the wave is then given by the typical wave number of the 
source. Note however, that at very high wavenumbers, 
(t&n — ti n )(k — LO s ) » 1 the last factor in Eq. (fil)|) al- 
ways leads to suppression. This behavior is illustrated in 
Fig.ffl 

Even though this flat space result cannot be taken over 
literally in a cosmological context, since we neglected the 
expansion of the universe, let us discuss in its lights, at 
least qualitatively, some examples of cosmological GW 
sources. 




FIG. 1: The quantity sin 2 ((fc - uo a )A)/(k - uj s ) 2 , A = (t fin - 
£in)/2, is shown as a function of k — uo a for the values of A = 
1 (black, solid) 0.3 (red dotted) and 10 (blue, dashed). For 
clarity the red (dotted) and black (solid) curves are displaced. 



If turbulence is generated in the early universe dur- 
ing a first order phase transition, as discussed in the in- 
troduction, one has the formation of a cascade of ed- 
dies. The largest eddies have a period comparable to 
the time duration of the turbulence itself (of the phase 
transition). They are not spatially correlated, and are 
described by a wavenumber power spectrum which is 
peaked at a scale corresponding to their size. Accord- 
ing to Eq. (fir?]) , these eddies generate GWs which inherit 
their wavenumber spectrum, since the last factor in this 
equation is flat enough, not to significantly influence the 
spectrum. Smaller eddies instead have higher frequen- 
cies, corresponding to periods smaller than the duration 
of the turbulence. According to Eq. (fTo]) , smaller ed- 
dies can imprint their frequency on the GW spectrum. 
However, they only contribute to the slope of the GW 
spectrum at high fc: the peak of the GW spectrum is 
always dominated by the largest eddies. Therefore, the 
GW spectrum is peaked at a wavenumber corresponding 
to the size of the largest eddies. 

Primordial magnetic fields are another cosmologi- 
cal source of GWs. They simply evolve by flux conserva- 
tion and dissipation of energy at small scales, and have no 
characteristic oscillation frequency. Their lifetime is typi- 
cally from their creation in the early universe until matter 
domination (after which only a negligible amount of GWs 
is produced [l(|)- Also in this case the wavenumber 
is imprinted in the induced GW spectrum. The suppres- 
sion for fc(ifm — tin) 3> 1 is seen by the fact that GWs are 
generated mainly on super- horizon scales, kt < 1. Once 
a scale enters the Hubble horizon, further GW genera- 
tion is negligible and the produced wave evolves simply 
by redshifting its frequency. 

A similar case is represented by GWs produced by 
the neutrino anisotropic stresses, which generate a 
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turbulent phase or simply evolve by neutrino free 
streaming [l^]. In the first case, the source is short lived 
and therefore its typical wavenumber is imprinted on the 
GW spectrum. In the second case there is no typical 
frequency and, like for magnetic fields, the relevant time 
is the time until a given wavenumber enters the horizon, 
kt ~ 1. Again, the GW frequency corresponds to the 
source's wave number. 

A counter example is the stochastic GW background 
produced by binaries of primordial black holes 
which are copiously produced in some models with extra- 
dimensions [131 ]. In this case, one has a well defined typ- 
ical frequency and the sources are long lived, so that 
w s (tfi n — ti n ) 3> 1. The GWs inherit the frequency of 
the source. 



IV. CONCLUSIONS 

To elucidate further our findings, let us consider again 
a spatially homogeneous source with well defined wave 
number k, which is alive only for a brief instant in time 
(a typical cosmological source). To determine the fre- 
quency of the generated GW, we observe it for some pe- 
riod of time (which is much longer than the inverse of 
the frequency). Since the source is active only for a very 
short period, as time goes on we observe GWs generated 
in different positions: what we see as happening in time, 
actually corresponds to a graviton-snapshot of the source 
in space, taken at the instant of time when it was on (see 
the bold solid lines in Fig. [5J). 

Conversely, for a source which occupies a very small 
spatial extension but lives for a long time (e.g. a binary 
star system), the situation is inverted. The GW signal 
arriving at the observer comes from virtually always the 
same position, but from different times. In this case the 
GW really monitors the time evolution of the source and 
therefore acquires its frequency (see the bold dashed lines 
in Fig. (2). 

This duality between k and uj is actually very natural, 
since they are related by Lorentz boosts. In general, both 
the spatial and time structure of a source lead to the 
generation of GWs. 

In this brief report we have shown that the standard 
lore that GWs oscillate with the same frequency as their 
source is not always correct. Special care is needed in the 
cosmological context, where the universe has the same 
temperature everywhere and can undergo some drastic 
change simultaneously (but in a uncorrelated way) in a 



brief interval of time. This can lead to sources of GWs 
which are very extended in space, but short lived. More- 
over, GWs from sources which only have a very slow time 
dependence and have no typical frequency, inherit usu- 
ally the spectrum of the source in fc-space. 

In general, a source can have a complicated spectrum 
with structure in both fc-space and frequency space. If 
the situation is not as clear cut as in the Gaussian wave 
packet, a detailed calculation is needed to obtain the cor- 
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FIG. 2: The past light cone (on which gravitons propagate) 
of a detector at x — observing for a given time interval 
sweeps a portion of space in the temporally confined source 
(solid, red line). Therefore, the wave number of the source is 
registered as the frequency of the emitted GW observed at x. 
On the other hand, the detector observes the time evolution of 
a spatially confined source (dashed, blue line), and therefore 
sees its frequency. 



rect GWs spectrum. On the top of that, as one sees ex- 
plicitely in Eq. ([5]), only the values (k, ui) with |k| = \u>\ 
of the Fourier transform of the source contribute to the 
GW amplitude. 
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